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Abstract 



Einstein's spherically symmetric interior gravitational equations are investigated. Fol- 
lowing Synge's procedure, the most general solution of the equations is furnished in 
case T\ and T\ are prescribed. The existence of a total mass function, M(r,t), 
is rigorously proved. Under suitable restrictions on the total mass function, the 
Schwarzschild mass M(r, t) = m, implicitly defines the boundary of the spherical 
body as r = B(t). Both Synge's junction conditions as well as the continuity of the 
second fundamental form are examined and solved in a general manner. The weak 
energy conditions for an arbitrary boost are also considered. The most general so- 
lution of the spherically symmetric anisotropic fluid model satisfying both junction 
conditions is furnished. In the final section, various exotic solutions are explored us- 
ing the developed scheme including gravitational instantons, interior T-domains and 
D-dimensional generalizations. 
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1 Introduction 

As motivation, let us consider various solutions of a toy model of the partial differential 
equation 

d 2 W(x,t) d 2 W(x,t) 



dx 2 dt 2 



To- (1) 



Here, To is a prescribed constant, (i) A particular solution is provided by W(x, t) = 
x — t + ^f{x 2 — t 2 ) which satisfies the initial value problem W (x, 0) = x + T x 2 /4 and 
9W gt'^ u_ Q = —1- (H) A class of general solutions of the same equation is given by W(x, t) = 
h(x + 1) + Tq/4(x 2 — t 2 ), where h is of class C 2 but otherwise arbitrary. This class contains 
infinitely many solutions but excludes infinitely many other solutions including the solution 
in (i). (in) The most general solution of this partial differential equation is furnished by 
W(x,t) = f(x - t) + g(x + t) + T /4(x 2 - t 2 ). Here both / and g are of class C 2 and 
otherwise arbitrary. This class contains all possible (smooth) solutions of the equation. 

Einstein's gravitational field equations, + 87r^|T*- = 0, inside matter are a 
system of second order, quasilinear, coupled partial differential equations in four space- 
time variables. It is almost impossible to obtain the most general solution for such a 
system in case the T*-'s are prescribed. However, if the space-time admits a group of 
motions or symmetry, then the equations simplify considerably. In fact, in the arena of 
spherical symmetry, using the curvature coordinates, Synge [TJ obtained the most general 
solutions where T\ and T\ are prescribed (the most logical prescription from a physical 
perspective). The interior was continuously patched to the exterior Schwarzschild metric 
across the junction of the spherical material in the local sense. However, the mathematical 
conditions assuring the existence of a boundary were not derived. Moreover, satisfaction 
of Synge's own junction conditions, T l ^ dD was not completed. 

In section 2, we write the spherically symmetric interior equations in curvature 
coordinates. Then, we exhibit the most general solution following Synge's prescriptions. 

In section 3, we prove the mathematical existence of a function, M(r, t). Physically, 
this is the "total mass" of the spherical body with coordinate radius r at coordinate time 
t. Under some reasonable assumptions, the implicit function theorem |2J guarantees the 
existence of a solution to r = B(t) for the equation M(r,t) = m, the Schwarzschild mass. 
The curve r = B(t) yields, in a natural way, the desired boundary for the spherical body. It 
is important to note that this patching is general and is therefore valid for junctions between 
various interior layers (as in, for example, multi-layered stars) as well as interior-vacuum 
patching. 

In section 4, we obtain necessary and sufficient conditions for the satisfaction of 
Synge's junction conditions pQ across the junction. Moreover, we also investigate the Israel- 
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Sen-Lanczos-Darmois (ISLD) junction conditions |3] across the junction and obtain general 
solutions of the problem. 

In the next section, we examine the weak energy conditions [3] thoroughly for the 
spherically symmetric scenario. We obtain the general solution of the inequalities in terms 
of four arbitrary slack functions. 

In section 6, the class of spherically symmetric [T*-] with real eigenvalues is criti- 
cally studied. As a particular application, the anisotropic fluid model (which contains the 
perfect fluid as a special case) is explored exhaustively. Theorems are proved on the most 
general solution of the corresponding field equations with both junction conditions of Synge 
and those of ISLD. Other special examples (black holes etc.) are also treated. 

In the last section, exotic spherically symmetric solutions and their relation to 
the proposed scheme are explored. Signature changing metrics as well as the Euclidean 
gravitational instantons j^j are furnished. Next, T-domain || equations and general so- 
lutions are provided. A special class of T-domain solutions yields the so called eternal 
black holes. Another special class of T-domain solutions involve complex eigenvalues of the 
stress-energy tensor. Such examples were already found in exotic black holes jjj. Finally, 
we give motivation for, and briefly investigate, spherically symmetric interior equations in 
arbitrary dimension D > 3. The corresponding general solution is provided jS]. 



2 Solution of the spherically symmetric field equa- 
tions 

We adopt notations and conventions from Synge's book [Tj, except that covariant 
derivatives are denoted by Physical units are chosen so that c = 1 and k := 8irG. 
Einstein's gravitational equations are furnished by: 

Sij '.— Gij ~\~ KTij = 0, (2a) 
T l := VjT ij = 0, (2b) 
VjS tj - kT { = 0. (2c) 

It is assumed that the metric functions, gij(x), are of class C 3 and the functions Ty(x) are 
of class C 1 . 

A spherically symmetric metric, in the curvature coordinate chart, and the natural 
orthonormal tetrad are characterized by 

ds 2 = e a ™ dr 2 + r 2 [d6 2 + sin 2 9 d^ 2 ] - dt 2 , (3) 
e\ l} = e\ 2) = r"^ 2) , ej 3) = (r sin 9)'%^ e\ A) = e^%. 
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Non-trivial equations and identities from (j2aHHJ) are provided by: 



g\ =r- 2 [1 - e~ a (l + 7-7,1)] + kT\ = 0, 



-7,ii + 2~ ( r 7,i + 2) (a - 7X1 



a 44 + -a,4 (a - 7) ,4 



c 4 



[1- 



7i =T 1 1>1 + T\ A 

1 2 

2 - ' /• 



re- a ),i] + 
2 



r 



r 



1 + 47,1 



T\ + -{a + 1 ) A T\ 



- ^l,iT\ - -Tl = 0, 



T 4 =T 4jl + T 4 4 )4 + 2 - [l + £(a + 7 ),i] ^4 + ^«,4 (T 4 4 - Ti) = 0, 
£i,i+£ 4 i,4 + 7 [l + ^7,i] *i + 5 (« + 7),4^1 



^7,1^4 _ = °' 



1 

2 ■ - r 

+ - [l + ~ A (a + 7),il £4 + ^a,4 (£ 4 4 - £1) - k% = 0. 



r l 4 v "'M * 2 
We study and solve these equations in a two dimensional domain given by: 

D := {(r,t) : < r < B(t), t l <t<t 2 \ . 



(4a) 

(4b) 

(4c) 
(4d) 

(4e) 
(4f) 

(4g) 
(4h) 

(5) 



Note that one may relax the restriction to the domain r < r < B(t) in which case radial 
integrals in the following should possess the lower limit of r . In such a case, an inner 
boundary will exist at r and the junction conditions discussed later should be applied 
to the inner boundary as well. The outer boundary curve, r = B(t), will be explicitly 
determined later (see fig. 1). 

Synge's strategy of solving the field equations is the following: 

• Prescribe T\ = from desirable physical properties and solve £\ = to obtain 
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Prescribe T\ = T }1 or relate it to T\ by an equation of state and solve £\ — £\ = 
to obtain e 7 = —gu- 

Define T\ by the equation £\ = 0. 

Define T\ = T^2 by the conservation equation T\ = 0. 
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By the preceding step, the identity Pg| ) implies that £ \ = 0. 

By the identity (I4h|). the conservation equation T4 = is satisfied. 



At this stage, all the field equations, conservation laws and identities are satisfied. One may 
impose further restrictions to the above scheme. For example, in the case of the perfect 
fluid, the conservation equation (jjej) becomes a differential equation for the pressure (or 
the energy density, if an equation of state exists), which must be solved. As well, one may 
require that further equations, such as matter field equations, need to be satisfied. 

Regardless of the variants on the above scheme, all solutions must satisfy the 
following most general solution yielded by: 



(6a) 

(6b) 

(6c) 
(6d) 



e -<*(r,t) 


=1 - 


- m- 

r |_ 




=e~ a 


( p '*> jexp 


T\(r,t) 


1 


Vit)f[t) 


Tl = Tl 





T 4 (x, t) x 2 dx 

rr 

h(t) + K / [T\(x,t) -T\(x,t)] e^'^xdx 



'0+ 



1 + 47,1 



T\ + \{a + i) A Tt- r -p tX T 



with f(t) :-- 



dm 

dt 



. Here, f(t) and h(t) are two arbitrary functions of integration which are 



of class C 3 . Synge lj set f 2 {t) = to avoid a singularity at the center. However, this 
function may be important in certain cases such as the study of wormholes. The function 
h(t) was absorbed by a transformation of the time coordinate though this is not always 
possible jH] - [HI- We retain these functions for generality and to satisfy junction conditions 
later. 



3 Conservation equations, the total mass function and 
the boundary 

We notice from the equations (jlc|) and (|4djl the existence of two additional differential 
identities: 



(r 2 G\) 1+ (r 2 G\)^0, 
{a + j) A G\ + a A (G\-G\) = 0. 



(7a) 
(7b) 
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However, because of V^G \ = 0, only one of the above additional identities is independent. 
Therefore, there must exist additional conservation equations 



%a :=(r 2 n) 1 + (r 2 r 4 4 ) i4 = 0, 

% b :={a + j) tl T\ + a 4 (T 4 4 - T\) = 0, 



d 



J r 2 T\{r, t) dt - ^r 2 Tl{r, t) dr 



0. 



(8a) 
(8b) 
(8c) 

(8d) 



The first of these equations has a simple physical interpretation. Integrating over a sphere, 
the equation relates the rate of change of energy in a sphere of radius r to the total energy 
flux entering or leaving the boundary of the sphere. The equation (jHaj) in the star-shaped 
domain D guarantees, by the converse Poincare lemma the existence of a function 
M(r, t) of class at least C 2 such that 



dM{r,t) 
Afi 



KT 



[T\{r,t) dt-T\{r,t) dr] 



\r 2 T\{r,t). 



From the equations (J4HJ), (pjaj) and (jHc|) . we conclude that 

2M(r, t) =/t/ 2 (t) - k / T^(x, t)x 2 dx, 

Jo+ 



hm M(r,t) =^/ 2 (t), 



-o(r,t) 



2M(r,t) 
r 

2M(r,t) 



exp [/i(t) 



X(r,t) :=/c 



x - 2M(x,t) 



X(r,t)} 
x 2 dx. 



(9a) 
(9b) 
(9c) 

(10a) 
(10b) 
(10c) 
(lOd) 

(lOe) 



We tacitly assume that r — 2M(r, t) ^ in D. The physical interpretation of M(r, t) is the 
"total mass" contained in the spherical volume of "radius" r and at "time" t. 

Next we wish to study the level curves of the function M(r,t). For the existence 
of such curves we state the following version of the implicit function theorem [2|. 

Theorem: 1. Let M(r,t) be a function of at least class C 1 in D such that for a point 
(r ,to) D, the function M(r ,t ) = c, a constant. Suppose that M jl |( r0jto ) ^ 0. Then 
there exists a function B(t; c) of class at least C 1 in the neighborhood of (ro,to) such that 
r = B(t; c) is a solution of M(r, t) = c in that neighborhood with r = B(to; c). 
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The boundary curve r = B(t) of the spherical body in the definition (J3J) is defined 
by the following: 



B{t) := lim B(t;c), 

c-^-m— 

3D := {(r,t) : r = B(t), h < t < t 2 } 



'IT 



(see figure 1). Here, m > physically represents the total Schwarzschild mass of the body. 
It is assumed that M > 0, Mi > or T\ < in D. 

t 

A 




boundary of D 



Figure 1: The considered domain with boundary, dD, which separates the interior domain D 
and the vacuum domain Dq. 



It is clear from the implicit definition of the boundary curve, dD : r = B(t) that 

M(B(t),t) = m, 
B(t) = 



~M A ~ 




r T i- 

1 4 




_Mi_ 


\8D 


rpA 
A 4. 


\8D 



(12) 



The spherical body collapses in case M i4 (r, i) > and expands in case M^{r,t) < 0. 

In case the measurable speed of the boundary is less than the speed of light, we 
must have: 



"\dD 



B(t) 



[e-^M^-e-^M^) 2 } 



dD 



rpl 

1 4 



< 1, 

> 0, 

< 



(13) 



\8D 



2m 



exp [h(t) + X (B(t),t)} 



B(t) + 2m. 



The interior domain D, the boundary dD and the exterior (vacuum) domain D 
are explained in the equations (Jo f a )l . (|6b"j) and (fTT]) and in the figure 1. Following Synge, we 
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shall now match continuously the interior metric to the exterior metric (transformable to 
the Schwarzschild chart). We must use the equations (|10d|) . (jlOeJ) and (fTTj) to arrive at: 



9 L1 (r,t) 



-ce(r,t) 



1 - 



1 



2M(r,t) 
r 

2m 
B(t) 
2m 



for < r < B(t), h<t< t 2 , 

for r = B{t), t\<t< t 2 , 

for B{t) < r < oo, h < t < t 2 , 



(14) 



j _ 2M(r,t) 
1 



exp [h(t) + x(r, t)] for < r < £ x < t < t 2 , 



2m 
B(t) 



exp [h(t) + X (B(t),t)] for r = B{t), t x <t <t 2 , 

1 - 2™] exp [&(t) + x(5(t),t)] for B(t) < r < oo, t 1 <t<t 2 . 

(15) 

The exterior metric can be easily transformed to the Schwarzschild coordinates. We shall 
next investigate both Synge's P and ISLD's [Hj junction conditions. 



4 Junction conditions 



4.1 Synge's junction condition 

Synge's junction conditions read 

T\n 



r H\dD 
TLin % 



0. 
1. 



(16) 



with n % a unit normal at the boundary. In the present case, with the help of (|13|) we can 
write for the relevant normal components: 



Hi 



Mi 



v/e-^Mi) 2 -e-T(M 4 ) 2 
M 4 

yV^Mx) 2 - e"^(M 4 ) 2 



3D- 



;i7) 



The equations (|Tfi|) reduce to 



[T 1 1 M 1 +T 4 1 M 4 ] |eD = 0, 
[T 1 4 M, 1 +r*M, 4 ] |8i3 = 0. 



(18a) 
(18b) 



By the equations (|9bjl and (|9cj). the junction condition ()18b|) is identically satisfied. More- 
over, the other junction condition (|18a|) yields: 



1 1 1 4 







dD 
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or, 



[e^ a T\M jl -T\M A ] ldD = 0. (19) 

There exist two possible cases here. In case the boundary is static, we must have from ()12j) 
and © 

B{t) = 0, 

T\Tt ldD = 0, (20) 
h{t) is an arbitrary function. (21) 

This case does not imply that the interior metric is necessarily static. 

In case the boundary is non-static, we obtain from (J9b"|) . (}9*cj) . (|T2^) and (fT9"|) 

B(t) ^ 0, 



IrpA 



T\T 



4\dD 



< 0, 



1 - 



2m 



-2 



exp [-*(£(*),*)] B{t) 



n 2 



t4 
1 4 



> 0. 



(22) 



\8D 



Thus, the function h(t), which originated as an arbitrary function of integration, can be 
utilized to satisfy the junction conditions. 

4.2 Israel-Sen-Lanczos-Darmois junction condition 

Next we consider the ISLD junction conditions. Namely, we consider the continuity 
of the second fundamental form at the junction. For this purpose, the three-dimensional 
metric for the hyper-surface corresponding to 3D x S 2 is obtained from (|15|) as 

B 2 (t) [oie 2 + sin 2 e d<f) 2 ] 



da 2 :— ds 2 r=B ^ 



1 - 



2m 

W) 



exj>[h(t)+ X (B(t),t)}- B{t) \dt 2 . (23) 



The extrinsic curvature components ^T] calculated from the interior and exterior metrics 
are the following: 



K, 



± 



lim 

(5^0+ 



re 



-a _ p -7 



B(t) 



\r=B(t)±S 



sin 2 9K± 



(24) 



2K t 



± 



lim 



B{t) 



2B{t) + e 7 - a 7>1 + B{t) {2a - 7) 



+ 



Bit) { a -2i) x - B{t) 



e a ~^a 4 



\r=B(t)±S 



It is clear from (jT4")l . (fT5|) and that 

To show the continuity of K u across dD, we consider the function 



(25) 
(26) 



-i 2 



1 3 



e a ~^a 4 



2L± := lim <^ e 7 - Q 7,i + 5(t)(2a - 7 ) i4 + £(t) (a - 2 7 ),i - £(t) 

I L J L J ) r=B(t)±S 

(27) 

The continuity of L across 3D implies, from (fTK|) and (|4al - l4fjl . (after a long calculation) 
the following algebraic equation: 







2m 



2k" 1 . 

B(t) 

U(t)e 2h ® + V(t)e h(t) + W(t). 



e h(t)+x{B(t),i) j^- _ £+j 



(2? 



Here, 



V(f) 
W(t) 



B(t) 



1 - 



2m 

W) 



1 2 



2 x (B(t),tW 



-B{t)B\t) { e x( fl (*).*) (Tl-lt)} 



|r=fl(t) 



(29) 



1 - 



2m 



T"l4 

J 4|r=B(t)- 



Analyzing the above quadratic (or possibly linear) equation for e h<yt \ we obtain the following 
solutions: 



Case I: T\, dD = B{t) T 4 , dD = and h(t) is arbitrary; 



Case II: e h{t) 



B(t) 



i 2 



-x(B(t),t) 



2m 

W) 



1 -2 



for 



(30a) 

i) T\ {dD = 0, B{t)T\ dD ± 0; 
it) (T\ - T\) 9D = 0; 
Hi) T\, dD ^ 0, Ti« n = 



Case III: T\ m > 0, T% dD < 0, B(t) ^ 0, 



L|(9Z) 

and e fc(t) 



B(t) 



2 



-x(B(*),t) 



1 - 



2m 



-2 



ry-i4 

1 4 



4|9D ~~ w > 

(30b) 



(30c) 



I9D 



It is clear that Synge's conditions ()21)1 and ()22)1 satisfy the ISLD conditions ()30aj) and 
IcJ). Case II represents a possible mathematical extension of the ISLD junction conditions 
to a non-time-like boundary. 
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5 Weak energy conditions 

Next the weak energy conditions in spherical symmetry are studied. We consider an 
observer with an arbitrary boost which, to our knowledge, has not been calculated before. 

In terms of the orthonormal components (denoted by indices in parentheses), the 
weak energy conditions jlj can be stated as 

T [am u^u^ > (31) 

for every time-like vector u'"' satisfying 

[ U W] 2 + 2 + [ u (3)] 2 _ [ u (4)] 2 = _ 1; (32) 

with > as dictated by reasonable physics. 

The general solution of the above non-linear algebraic equation ()32j) is given by: 

= sinh/3cos#, vP^ = sinh/3sin#cos0, = sinh (3 sin 9 sin 0, w^ 4 -* = cosh/3, 
P e R, 9 e (0, tt), e (-tt, tt). (33) 

For spherical symmetry, choosing the orthonormal basis of Q, the inequality (J31|) together 
with equations (JBlIj) yield: 

[T ( i)(i) - T (2 )(2)] sV + ^(2)(2)a; 2 + 2T (1)(4 )xy + T (4 )( 4 ) > 0, (34) 
x:=tanh/9, y := cos9. 

Analyzing the inequality for all (x, y) G [—1,1] x [—1,1], we conclude, after much 
calculation, that either 

(i) T{i)(i) = ^(2)(2), ^(i)(4) = 0, T(4)( 4 ) > 0, T(4)( 4 ) + T( 1)(1 ) > 0, (35a) 
or else 

(ii) T (1)(1) > T (2)(2) , (T (1)(4) ) 2 < T (4)(4) [T (1){1) - T (2)(2) ] , 

P(l)(4)) < [T(l)(l) -T(2)(2)] [T( 4 )(4) +^(2)(2)] • (35b) 

We can solve the inequalities (|35aj) ()35b|) by utilizing slack functions: 

T\=T§ = -E 2 (r,t), (36a) 

T\ =T% = P 2 (r,t) - [E(r,t)smQ(r,t)} 2 , (36b) 

T\ =Tg = P 2 (r,i) - [£(r, t) sin Q(r, t)] 2 + tf 2 (r,t), (36c) 
e («-7)/2 T i =T(i)(4) = _ jff(r; t)jE(r? t) cos Q(r, t) (36d) 

Here, the slack functions E(r,t), P(r,t), Q(r,t), H(r,t) are of class C 1 but otherwise 
arbitrary. 
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6 Real eigenvalues of T l - and anisotropic fluid mod- 
els 



First, we analyze and solve the problem of a spherically symmetric T* possessing real 
eigenvalues. Recall that the eigenvalue problem for T* is given by: 



T ) E U = X ^ E U- 



(37) 



In the spherically symmetric case, T\ = T\ = T\ = T\ = 0, T\ = T\. Therefore, the 
eigenvalues of T % - are given by: 



A(2) 


= A( 3 ) = Tl, 




(3 


2A ( i) 


= T\ + T\ + \ 


/A, 




2 A ( 4 ) 


= T 1 1 + T 4 4 -^ 


/A, 




A : 


= {t\-t\) 2 


-e Q " 7 (T\) 2 . 





It is clear that A < will imply complex eigenvalues. We restrict ourselves in this section 
to the case where the stress-energy tensor possesses real eigenvalues (A > 0). 

In a static model, T\ = and A > is automatically valid. In case of the weak 
energy conditions in (|35a)L (|35b|) and the corresponding solutions in ()36al - IBlxIj) . 

(39) 



A = P A + 2P 2 (H 2 + E 2 cos 2 Q) + (H- E 2 cos 2 Qf > 0, 



and thus real eigenvalues are guaranteed. (However, A > may not imply the weak energy 
conditions.) 

Assuming the existence of real eigenvalues, the corresponding natural orthonormal 
eigenvectors are furnished by: 



E h) 
Eft) 

E h 



r °(2)'- C/ (3) 



(rsinfl) E} 1} = 



E (i) 



o, 4) 



'(3): ^(1) 
'(4) 



V{1) T 1> E (A) = -V(4)T\, E {A) 



T\ — T\ — VA 

2 _ Z73 — Q 



E k 



(40) 



hi)]' 



2e" 



y/A T\ - T\ + y/A 

2e~ 7 



y/A T\-T\ + VA 



The decomposition of the spherically symmetric T % i in terms of the real eigenvalues 
and eigenvectors is accomplished by 



[A(i) - A (2) ] E l {1) E J (1) + \ {2 )9 V + [A( 2 ) - A (4) ] E l {A) E J {4) . 
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(41) 



It is worth noting that the above algebraic structure of the stress-energy tensor is common 
to many different physical arenas. For example, the anisotropic fluid is specified by 



Pn '■= P± := A( 2 ) = A (3) = T 2 , n := -A( 4 ), 

u % := -E(4), s l := -E(i), 
T ij = (v + p±)u i u j +p ± g ij + (p l -p ± )s i si. (42) 

The physical quantities are the energy density (/x) the radial pressure (p„) and the angular 
or transverse pressures (p±)- Anisotropic fluid models have received much attention mainly 
in the arenas of stellar structure theory, black holes and cosmology [7j, [T2j, [131, C3- Note 
that the nomenclature "anisotropic fluid" is misleading. The stress-energy tensor in ()42j) 
actually represents a fluid which is not necessarily isotropic. 

In case Am = A( 2 ) = A( 3 ) (or p n = p±_ =: p), the equation (j%2j) yields the well known 
perfect fluid stress-energy tensor: 

T ij = (fx + p) uV + pg ij . (43) 

This equation implies, by equation (|65)l and ()38|). the isotropy equation: 

9 

T* I- - 7" 

- + ^ (Ti - T 4 4 ) 



(Ti - T 4 4 ) 2 - e 2a - h -*(T\) 2 = rT l xl 



(T\ - It) 

_ e -(^+x)/2 j- e 2a-(h+x)/2 T r| (44) 



.4 

It is a formidable equation to solve in general (see |T3] for detailed considerations of the 
static case.) 

In case the spatial eigenvalues are identically zero, the stress-energy tensor in (J42|) 
reduces to that of an incoherent dust. 

In case we identify p = A (2) = A(3), /J = — A(4), a := Am — A (2), the stress-energy 
tensor is, 

T ij = (/i + p)«V + pg ij + as l s j . (45) 

The above T u is due to a perfect fluid plus a tachyonic (space-like) dust. Such a stress- 
energy tensor has been considered in a cosmological model ^31 where the dust contributes 
to the dark matter or dark energy component of the universe. 

Now we are in a position to state and prove the main theorems of this section 
involving anisotropic fluids. 



Theorem: 2. Let the spherically symmetric interior equations J^aH^o!| ) and the conserva- 
tion equations fije\ ), ( pjty hold in the coordinate convex domain D defined by ([3J). Moreover, 
let the stress-energy tensor T*- be that of an anisotropic fluid given by Also, let the 

physical conditions T\ < and T\ — T\ > be satisfied in D. Then, the most general 
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solutions of all the equations and inequalities are furnished by the following: 



< gi < 1 


, < q 2 < 1, 




2M(r,t) 


:= nqi<yF 2 {q 2 t) + J E 2 (x, q 2 t)x 2 dxj > 0, 




2M(r,t) 
r 




e T(n*) 


= e-^exp [/*(*) + X (r,t)], 




X(r,t) 


:= « ft / e*^ [£ 2 (x,g 2 t)cos 2 Q(a;, 


q 2 t) + P 2 (x,q 2 t)] xdx, 


A(r,t) 


- ( gi ) 2 ^ 2cOS ^ + p2 ] 2 2 4 4 eQ " 


" 7 (M i4 ) 2 > 0, 


—fi{r,t) 


= (g 1 )- 1 v / A + ,E 2 [l + sin 2 Q] - P 2 


>o, 


2 , x 
—Pu{r,t) 


= (q 1 )- 1 VX-E 2 [l + sin 2 Q] +P 2 


>o, 



(46) 



= ^ [r2( p 2 _^ sin2g)]i 

+r | i(E 2 cos 2 g + P 2 ) 7 ,i - (/c<?i) _1 e~ (a+7)/2 [ e ( 3Q -^ 2 M, 4 ] A 

Kp7 /2 r 2 Al/4 I — ' " U - U > 



u 4 



s 1 



^^1/4 cos2 Q + p2 ) + ^ > °. ( 4? ) 

V^ 2 cos 2 Q + P 2 )^, s 2 = 0, 



. 2V2e a ^M 4 
s = =F 



Kr 2 A l/4^ gi (£2 cog 2 Q + p 2 ) _ 

ifere, the functions F(q 2 t), h(t), E(r,q 2 t) (not identically zero) are of at least class C 3 in 
D. Aside from these restrictions, the functions are arbitrary. 

For proof of the above theorem we used the equations (j6aH6d|) . (jlOal - ITUej) . (jHfiai - 
(jHEJ), Ptljl . (jUj) and P^jl . The two parameters in (jUjj) may appear to be superfluous. 



However, note that in the limit — > 0+, the solutions in (j4*U|) yield the flat space metric. 
Moreover, for < q\ < 1 and and limg 2 — * 0+, the metric goes over to a static one. 
Furthermore, sufficiently small positive values of q\ and q 2 facilitate satisfaction of the 
complicated inequalities 1 — — > and A > 0. 



14 



We consider here a specific example of an exotic black hole. Consider the following: 



F(q 2 t) 

2M(r,t) 
e -a(r,t) 

T\(r,q 2 t) 
X(r, q 2 t) 



= 0, E\r,q 2 t) : 



1 



l-q 2 t 



3< j. 

(2m 
\KQlJ 



1 - q 2 ty ' 

j /2m\ 1/j 



j'-i 



k- 2 qi E 2 (r,q 2 t), V3<k, 



(48) 



1 + A;" 2 ) In 



nqir 



i 



Kq±r 



s fc(*) 



The above describes a mathematically rigorous collapse model for an anisotropic fluid black 
hole 0. 

Now we shall consider the junction conditions for the solutions given in (|48p . We 
state and prove the following corollary to the preceding theorem. 

Corollary:. Let the conditions stated in the previous theorem be valid in D with E(r, q 2 t) ^ 
0. Moreover, let both Synge's junction conditions T^Mj\gu = and the ISLD junctions 
conditions [Kij]g D = hold on dD. Then, either, 

sin 2 Q(r, q 2 t) E 2 {r, q 2 t) =P 2 {r, q 2 t) + [B{t) - rf N{r, q 2 t) > 0, v 2 G {1, 2} U [3, oo) 



and 



rt 2 f rB(t) q ] 
F 2 (q 2 t) =(c ) 2 + J drlj x 2 — [E 2 (x, q 2 r)} dx\>0, 



or else 



/■*a f fB(t) o 

F 2 (q 2 t) =f(q 2 t) + jf * ^ x 2 ^ [£ 2 (x, g 2 r)] rfx 



(49a) 



(49b) 



£ 2 (r,g 2 t) 



P 2 



sec/i [P(r, g 2 i)] + sin Q 



> 0, f{q 2 t) ^ 0, 



and 



e Ht) =4 J n _ 



2m, 



-2 e - X (r,t) 



M 4 cosh(P(r, g 2 i)) 
Kqir 2 E 2 (r, q 2 t) 



> 0. 



|r=B(i) 



Pere, N(r,q 2 t), R(r,q 2 t) and f(q 2 t) are functions of at least class C 3 in D but otherwise 
arbitrary. 
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Proof: By the equations ()36c|) and (|49ajl it follows that 



[B{t)-r] v N(r,q 2 t), T 



■l 



= 0. 



l|r=B(t) 



Furthermore, 



2M 4 



dF 2 



(?2*) 



B(t) 



+ 



x 2 —E 2 (x, q 2 t) dx = 0. 



. K <?1 J |r=B(t) 



0+ 



Therefore, by (J2T)) and (fHOaj) . both Synge's condition and the ISLD conditions are satisfied. 
In the second case, by the equations in f)36c|) and ()49b|) it is deduced that 



K- l T\(r, q 2 t) = P 2 - E 2 (r, q 2 t) sin 2 Q = E 2 (r, q 2 t)sech 2 [R(r, q 2 t)} > 0. (50) 
Moreover, M 4 ^ and 



Thus, both equations (j22j) and ()30c|) are satisfied. ■ 

We have previously proved jS] that under the two conditions T\ = and 4 = 0, 
the solutions of the equations (|4aH4tj) can be transformed into a static solution. This is 
the interior version Birkhoff's theorem. We can investigate directly the static limit of 
equations ()4aH4tj) . Under suitable assumptions, including > q, the boundary, dD, of 

the spherical body is given by r = b, a positive constant. Now, the general solution will be 
furnished in the following statement: 

Theorem: 3. Let the static version of the spherically symmetric field equations and one 
conservation law ^4a\j4J\ ) hold in the domain D := {(r, t) : < r < b, t± <t < t 2 }. More- 
over, let the stress-energy tensor be given by fcfify , satisfying fi{r) > 0, /i(r) + p u > 0. If, 
in addition, both Synge's and the ISLD junction conditions hold at r = b, then the general 
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solutions of the static equations are furnished by: 

< qi < 1 , b > 0, c e R, v 2 e {1,2} U [3,oo), 



2M(r) 

A(r) 
P± 

(sin 2 Q(r))P 2 (r) 



/v 



(c ) 2 + / x 2 jx{x) dx 



> 0, 



1 - 



2M(r) 



k / e a ^ +Pn( a; )] xdx, 
Jo+ 

e x{r)-a{r) 

Hr)+P„(r)} 2 >0, 

giE 2 (r) > 0, p, = q 1 [P 2 (r) - (sin 2 Q(r)) E 2 



(51) 



— V Pn.l + ^"7,1 



2r 

e -7(r)/2 5 



(4)' 



±e a(r)/2« 



(1)' 



P 2 (r) + (b-r) v N(r) > 0. 



i/ere ; E(r), P(r), Q(r) and N(r) are functions of at least class C 3 . Moreover, these 
functions and the parameters cq and v are arbitrary save for the restrictions imposed above. 
Proof follows from (|47)|) and (|49ajl . Note that to avoid a singularity at r = (if this is 
included in the domain), the constant Co should be set equal to zero so that lim r ^ M(r) = 0. 
An illustrative example will be provided in the following: 



< qi < 1, 



E 2 (r) 



3. 



3 sin 2 Q(r) 



b > 0, c = 0, 
2M(r) = /«gir 3 , 
P 2 (r) - 3nqi(b - r) x 



6 + r 



Hir) 
p(b) 



3y/l — nqib 2 — a/1 — nqir 2 
1 - nq^ 2 , 
3yl — fi^fr 2 — \J\ — K^r 2 



1 — /tgir 2 + a/1 — /«gi6 2 



> 0, 



3^1 - Kgi6 2 - 1 

3<?i, p„ = P_l =: p(r) = 3gi 
0. 



1 — Kgxr 2 



1 — ftq^fe 2 



3>/l — Kgi6 2 — yT— fi^i? 



(52) 



The above obviously yields the well known interior Schwarzschild constant density solution. 

One final example which illustrates the use of this scheme is that of the inner layers 
of a static neutron star [TH]. In this case, the energy density, /x = 7(4) (4) = 



T4 is known 
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from the quantum mechanics of degenerate Fermions. As well, there is the ultra-relativistic 
fluid equation of state, which should be valid in the inner layers of the star. We summarize 
as follows: 



fi(ki 



S7T 



kp 



k 2 (k 2 + m 2 n ) 1 ' 2 dk 



Here h is Planck's constant, kp is the Fermi- momentum and m n is the neutron mass. Since 
the extreme relativistic limit is employed, the mass terms may be neglected compared to 
the Fermi momentum so that a pressure calculation gives: 



Vn = P± 



(53) 



Also, fj,(r) is obtained by utilizing along with the linear combination G\— G\ = |«/^(r). 
Now the isotropy equation of the general solution reads, in terms of /i(r), 



' 2M(r)~ 


-l - 


r 





+ M(r) 



(54) 



Noting relation (J9bj) and assuming a series solution in r we arrive at the following: 



M(r) 

which, from (jUaf) and (|6b|) yields: 
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r, /i(r) 



4 r 
7r 



3 

7kx 2 ' 

7 

4 



(55) 



with r a constant. The neighborhood about r = is excised as the singularity at this 
point is due to the ultra-relativistic approximation. Also, this solution is not valid for outer 
layers of the star where deviations from the ultra-relativistic case are significant. 



7 Exotic spherically symmetric solutions 



The exact solutions in (|10al - llOef) can be generalized by abandoning the weak energy 
conditions to express: 



rr 

2K- l M(r,t) =f 2 (t) - / Tl(x,t)x 2 dx, 

Jo+ 



e l{r,t) _ 


1 - 


rr 

X{r,t) :=k \ 






H(t) ± 0. 





r 

2M(r,t) 



e*%(t), 



(56a) 

(56b) 

(56c) 

(56d) 
(56e) 
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There are many situations when the solutions to these equations may prove to 
be "exotic" in some sense. For example, the equations (|56bH56ej) reveal that the condi- 
tion signature [g^] = +2 may not be preserved everywhere. A simple example may be 
considered in the exact vacuum metric given by: 



ds 2 



(l - ^ 1 dr 2 + r 2 (d6 2 + sin 2 9 d<p 2 ) + (l - ^ t 3 dt 2 



(57) 



In case t < 0, the metric is obviously transformable to the Schwarzschild solution. However, 
for t > the line element (J57|) yields the spherically symmetric vacuum gravitational 
instanton solution. In (JSTj) . lim^ 9u{ r , t) = 0, indicating the existence of a horizon. All 
the null rays from the Schwarzschild universe suddenly halt on such a horizon. It may be 
called the instanton horizon. Signature changing metrics in general relativity have been 
studied in [7j, ^7j. The most general spherically symmetric instanton solution in curvature 
coordinates is furnished by the equations (j5tial - l57)e|) with the choice H{t) = —e h ^ < 0. 

Next, consider spherically symmetric T-domain solutions ([6 and references therein). 
The metric is locally expressible as: 



ds 2 =e A ™ dR 2 + T 2 (d6 2 + sin 2 9 d<j) 2 ) - e v ™ dT 2 - 
D T := {(T,R): 7\ < T < T 2 , R 1 < R < R 2 } . 



(5? 



Einstein's field equations (j 1 +k;0*- = can be solved with the metric (|58|) . The general 
solution ("dual" to the solutions in (joal - l6djl ). is furnished by: 



1 

T 



a{R)-K / {T') 2 Q\{T',R) dT' 



To 



2~(T,R) 
T 



1 > 



(59a) 



2E(T,R) 
T 



exp ^(R) + k j T e" {T '> R) [9 4 4 (T', R) - 6\(T', R)] T' dT' 



e 2 = e 3 3 



1 + ? A, 4 



9 4 4 + j(A + ^) il © 4 1 -|a,4©\, 



(59b) 
(59c) 

(59d) 



all other 9Vs = 0. 



Here, the functions c(R) and (3{R) are of class C 3 but otherwise arbitrary. The "total 
tension" function, S(T, R) is generated by the tension density since, in the T-domain, it is 
@\ which appears in (|59aj) . This class of solutions includes eternal black hole solutions. 

Another special case of T-domain solutions occurs whenever the stress-energy ten- 
sor matrix [0 '•] admits complex eigenvalues. The algebraic criterion of such occurrence is 
provided by the strict inequality: 



a» := ( 4_e 1 1 ) 2 + 4e 1 4 e 4 1 <o. 



(60) 
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As an example, the following l - has appeared in the late stages of gravitational collapse 
studies [7]: 

< q 1 < 1, < q 2 < 1, 3 < j, v 7 ^ < k, 



e\(T,R)=- 



Q\(T,R) =k- 2 3gi 



1 - qzR) 3 



- < 0, 



1 - q 2 RV 



> 0, 



(61) 



q\(t,r) = mq2TJ ~ 2 . >o. 



Finally, we shall consider the spherically symmetric field equations in an arbitrary 
.D-dimensional manifold (with D > 3) [S|. There has been much study on the possibil- 
ity of extra dimensions in light of superstring theories. In the low energy sector, many 
of these theories reproduce a higher dimensional general relativity in which, above some 
energy scale, all dimensions may be considered non-compact. These higher dimensional 
field equations may therefore have relevance in these theories. 

The metric in curvature coordinates is provided by: 



Mr,t) dr 2 



r 2 dQ 2 D _ 2 ) 



5* r '*> dt 2 



(62) 



with 



D-3 



n=l 



\m=l 



D:={(r, 9 {0) , 
0<6 {0 ), 



■ , d(D-s), t) E R D : t x < t < t 2 , < n < r < r 2 , 
■,0(D-4) < 7T, < 6 {D _ 3) < 2tt} 



(63) 
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The D-dimensional field equations and conservation laws read: 
D-2 



2r 2 



[(D - 3) (1 - e- a ) - re- a 7i i] + kT\ = , 



£\ = — [7,401,4 - (a A f - 2a 44 ] 



, x2 2(£>-3), 
27,ii + (7,i) 2 + 1 - > (7 " «)l 



- 7 ,ia 1 + -( J D-3)( J D-4) 



+ 2p-3)P-4) + _ 



£ 4 



2r 

, -^3^ [P " 3) (1 - e- a ) + re- a a J + kT 



2r 2 

T 1 =T\ tl + T\a + 



1 = 



1 D-2 

o7,i + 

2 r 



2^ + r 



^ 4 =n 4 + n,i + pi - t d + \ t \ 



(a + 7)1 + 



2(D-2) 



= 0, 



(64a) 

(64b) 

(64c) 
(64d) 

0, 

(64e) 
(64f) 



The general solution of the Einstein field equations and conservation equations furnished 
utilizing the scheme in this paper is: 



e -a(r,t) =1 + 



2k 



(D - 2y 



D _ 3 j i 4 yx,t)x ox ^ D _ 3 -.1 ^ D _ 3 , 



7(r,t) 



-a{r,t) 



exp < h(t) + 



2k 



D-2 



n 



x fl-3_2M(x,t) 



(65a) 
(65b) 



Again, the functions /(£) and h(t) are of class C 3 but otherwise arbitrary. 



8 Concluding remarks 

In summary, the general solution to the spherically symmetric Einstein field equations 
was provided in the case when the energy density and parallel pressure are known. Both 
Synge's junction conditions and the Israel-Sen-Lancsoz-Darmois junction conditions have 
been studied and solved in general. The junction or boundary is defined by the existence 
of a total interior mass which has been rigorously proved in section 3. The weak energy 
conditions in spherical symmetry for arbitrary boost were presented and solved utilizing 
slack function methods. Specific matter models have also been considered including the 
anisotropic fluid satisfying both junction conditions, which includes the perfect fluid as a 
special case. Finally, exotic extensions were considered. 
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